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1 .  INTRODUCTION 


The  study  of  atmospheric  diffusion  can  be  simply  defined  as 
the  investigation  of  the  diluting  effects  of  the  atmosphere  on  a 
released  contaminant.  The  atmosphere  is  most  often  in  a 
turbulent  state  implying  that  its  motions  can  be  predicted  only 
in  a  stochastic  sense.  Since  virtually  all  gaseous  pollutants 
and  most  airborne  particulates  very  quickly  become  well  mixed  on 
the  molecular  level  in  the  troposphere,  diffusion  of  those 
materials  can  be  described  in  a  similar  probaba  1  i st i c  manner. 
Introducing  statistics  automatically  requires  the  user  to 
consider  the  variations  of  those  statistics  in  time  and  space. 

As  an  example  of  the  later,  a  predicted  crosswind  concentration 
distribution  represents  only  a  first  order  approximation  to  the 
actual  distribution,  which  will  often  show  significant  variations 
from  that  prediction.  The  study  of  those  fluctuations  is 
presently  only  beginning  (see  Sawford,  1985). 

This  study  focuses  on  the  time  variations  of  diffusion 
statistics,  specifically,  the  standard  deviation  of  the  crosswind 
concentration  distribution  (sigma-y).  Skupniewicz  and  Schacher 
(1984b)  have  shown  that  overwater  releases  of  material  will 
diffuse  in  the  crosswind  (y)  direction  with  two  scales  of  motion. 
On  one  scale,  material  will  disperse  about  the  centerline  of  the 
plume  or  cloud  due  to  turbulence  of  length  scales  close  to  the 
plume  or  cloud  size.  This  scale  of  dispersion  is  referred  to  as 
relative  dispersion,  implying  that  the  statistics  can  be  observed 


by  the  motions  of  parcels  relative  to  each  other  in  a  moving 
frame  of  reference.  On  the  second  and  much  larger  scale,  the 
plume  or  puff  will  disperse  relative  to  a  fixed  axis,  usually 
chosen  to  be  the  mean  wind  direction.  Turbulence  acts  on  the 
instantaneous  plume  or  puff  as  if  it  were  a  single  entity,  and  is 
referred  to  as  single-particle  diffusion.  A  more  common  and 
understandable  synonym  is  meander. 

Skupniewicz  and  Schacher  (1984a)  implicitly  considered  the 
combined  diffusion  properties  of  both  components  when  hourly 
averaged  plume  parameterizations  were  derived  from  tracer 
experiments.  In  this  report,  the  two  components  are  separated  in 
order  to  supply  more  information  on  the  nature  of  the  diffusion 
processes,  and  resultantly  more  knowledge  of  potential  hazard 
from  a  contaminant  release. 
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2.  EXPANDED  DATA  BASE 


Most  descriptions  of  plume  parameters  rely  on  empirical 
formulae  and/or  major  assumptions  about  the  physical  behavior  of 
the  atmosphere.  Verification  of  these  parameterizations 
therefore  should  utilize  data  gathered  under  a  wide  variety  of 
atmospheric  conditions  and  geographical  locations.  The 
Environmental  Physics  Group  (EPG)  at  the  Naval  Postgraduate 
School  (NPS)  has  made  a  significant  effort  to  collect  a  diverse, 
overwater  data  base  over  a  several  year  period  for  this  purpose. 
These  data  are  reported  in  Skupniewicz  and  Schacher  (1984a)  and 
will  be  hereafter  referred  to  as  the  NPS  data.  The  NPS  data 
primarily  consists  of  concentration  profiles  collected  from  an 
aircraft  platform  at  various  distances  from  a  continuous  surface 
release  of  SF6  gas.  Part  of  the  work  reported  here  was  to 
integrate  a  new  data  set  collected  by  the  German  Military 
Geophysical  Office  (GMGO)  (see  Groll  et  al,  1984)  into  our  data 
base.  These  data  were  obtained  in  the  North  Sea  via  continuous 
SF6  releases  from  a  ship  and  subsequent  downwind  plume  transects 
with  a  second  ship  sampling  gas  concentrations.  Table  1  lists 
information  on  the  three  experiments  made  available  to  EPG. 


Experiment  Data  No.  of  Sampling  Boats 


No.  of  Plume  Transects 


23 

May-06 

Sep 

79 

3 

516 

1  4 

Apr-29 

Apr 

80 

4 

558 

04 

Nov-1  3 

Nov 

80 

4 

260 

Table  1.  GMGO  overwater  tracer  experiments  available  at  EPG-NPS. 
Each  transect  produces  an  "instantaneous”  sigma-y  value  used  in 
this  analysis. 

Task  I  specifically  involved  several  steps  needed  to  produce 
plume  parameters  from  the  GMGO  data: 

1)  transferring  original  9-track  data  files  to  mass  storage 
files  in  IBM  3033  readable  format, 

2)  combining  half-hour  average  meteorological  files  with 
nearly  coincident  tracer  profiles,  and 

3)  calculating  the  second  moment  of  the  concentration 


distribution  (sigma-y)  for  each  profile  and  the 
downwind  distance  from  navigational  information. 
Table  2  lists  the  resulting  data  set  contents  and  format. 
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VARIABLE 

DEFINITION 

UNITS 

FORMAT 

DATE 

day,  month,  year 

no  units 

312 

TIME 

hour 

Z 

14 

RH 

relative  humidity 

% 

F5.1 

WD 

wind  direction 

deg 

14 

DB1  7 

dry  bulb  at  17m 

C 

F5.1 

WSM 

wind  speed  (MET) 

m/s 

F5.2 

WB1  7 

wet  bulb  at  17m 

C 

F5.1 

ET 

water  temperature 

C 

F5.1 

DB3 

dry  bulb  at  3m 

C 

F5. 1 

WB3 

wet  bulb  at  3m 

C 

F5.1 

GRAD 

gradient 

C/1  000m 

14 

SB 

speed  of  boat 

m/s 

F5.2 

SS 

source  strength 

rrn  /hr 

F2.0 

SH 

sampling  height 

m 

12 

WSS 

wind  speed  (SF6) 

m/s 

F5.2 

DS 

distance  from  source 

m 

F8 . 1 

BTS 

bearing  toward  source 

deg 

F5.0 

HD 

heading 

deg 

F2.0 

#  pt 

no.  of  points  used  for  CWCI 

no  units 

12 

CWCI 

cross  wind  concentration 
integration 

ppb-m 

F8 . 0 

MEAN 

mean  mass  position 

m 

F6 . 1 

SD 

standard  deviation 

m 

F6 . 1 

Table  2.  Contents  of  new  GMGO  overwater  tracer  data  set  added  to 
the  EPG  overwater  data  base.  Format  is  in  FORTRAN  code.  Public 
mass  storage  data  set  name  is  "MSS. F3896. DATA. NEWSM". 


3.  RELATIVE  DIFFUSION  PARAMETERIZATION 


The  basic  approach  used  in  parameterizing  the  relative 
diffusion  is  to  group  observations  according  to  Pasquill-Gifford 
equivalent  stability  classes  (see  Turner,  1967),  and  then  regress 
sigma-y  versus  range  (from  the  release  point)  for  each  class. 
Sigma-y  was  defined  in  the  last  section  as  the  second  moment  of 
the  concentration  distribution  for  a  given  profile.  The 
customary  assumption  that  concentration  in  the  y  direction  is 
independent  of  concentration  in  either  the  x  or  z  direction  leads 
to  the  conclusion  that  these  instantaneous  "snapshots"  of  the  SF6 
plume  are  identical  to  the  crosswind  dimensions  of  a  hypothetical 
"puff"  or  "burst"  release  under  the  same  atmospheric  conditions. 

The  methodology  for  determining  stability  class  is  described 
in  detail  in  Skupniewicz  and  Schacher  (1984a).  Briefly,  the 
scheme  uses  windspeed,  air-sea  temperature  differences  and 
relative  humidity  to  produce  stability  classes  equivalent  to  the 
well-known  Pasquill-Gifford  diffusion  categories.  [Whether  a  true 
equivalence  exists  is  a  matter  of  controversy  and  is  discussed  in 
Skupniewicz  and  Schacher  (1984a).] 

Under  this  scheme  classes  are  assigned  a  letter  designating 
atmospheric  stability,  with  "B"  representing  the  most  unstable 
situations,  "C"  representing  moderately  unstable,  "D"  representing 
neutral,  and  "E"  representing  moderately  stable  conditions. 

After  dividing  the  data  according  to  stability,  the  subsets 
were  examined  for  dependence  upon  wind  speed,  time  of  day,  and 
day  to  day  trends.  No  obvious  dependences  were  found.  Figures 


1-3  show  sigma-y  vs.  downwind  distance,  with  the  integer 
truncation  of  the  wind  speed  used  as  the  point  designator. 

In  as  much  as  windspeed  gives  a  first  order  approximation  to 
sea-state,  increased  surface  stress  within  a  given  stability 
class  has  a  negligible  effect  on  diffusion.  This  result  is 
particularly  surprising  for  the  neutral  case  (class  D),  where  a 
small  air-sea  temperature  difference  will  always  produce  neutral 
stability  under  this  scheme  regardless  of  windspeed.  Intuitively, 
one  would  expect  that  higher  windspeed  would  enhance  the  plume 
spread.  Because  stable  or  unstable  profiles  tend  towards 
neutrality  as  windspeed  increases  in  this  stability  scheme, 
windspeed  dependence  is  not  expected  in  non-neutral  categories. 

It  is  also  noteworthy  that  stable  situations  are  found 
exclusively  in  light  wind  situations  in  the  Worth  Sea  experiments. 
These  situations  were  not  found  in  the  California  experiments. 

Based  on  the  above  data,  the  parameterization  chosen  for 
relative  diffusion  is  a  simple  linear  relationship.  The 
following  factors  influenced  this  choice.  Instantaneous  releases 
into  a  turbulent  atmosphere  will  theoretically  experience  the 
following  four  different  growth  regimes: 

a)  the  initial  stage 

b)  the  inertial  stage 

c)  the  central  stage 

d)  the  final  stage 
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Figure  3.  Same  as  figure  1,  except  a)  Pasquill-Gif ford  stability 
class  D  (neutral),  NPS  data.  13  data  points  hidden,  145  data  points 
missing,  b)  class  E  (slightly  stable),  NPS  data.  1  data  point 
hidden,  45  data  points  missing. 
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The  initial  stage  applies  to  distances  of  only  10-IOOm.  The 
inertial  stage  will  be  relevent  from  100m  to  less  than  1km.  The 
central  stage  will  apply  from  1km  to  10km  and  the  final  stage 
will  apply  thereafter.  These  bounds  are  highly  variable. 

The  initial  stage  of  growth  is  valid  when  the  instantaneous 
puff  "remembers"  its  initial  size.  In  this  regime  cloud  size 
will  increase  linearly  with  distance  or  travel  time  because  each 
particle  in  the  cloud  will  move  with  its  initial  velocity  and  the 
particles  have  not  had  time  enough  to  become  correlated  with  each 
other . 

The  initial  stage  is  followed  by  growth  dominated  by 
turbulence  in  the  inertial  subrange.  In  the  inertial  region, 
plume  growth  can  be  shown  to  follow 

o2  =  aet^  ,  ( 1  ) 

yr 

where  oyr  is  the  relative  si gma-y , 
a  is  a  constant, 

e  is  the  turbulent  kinetic  energy  dissipation  rate,  and 
t  is  travel  time  of  the  cloud. 

This  inertial  stage  of  growth  is  valid  only  where  the  width  of 
the  cloud  is  small  compared  to  the  average  height  of  the  cloud; 
surface-released  clouds  seldom  meet  this  criterium. 


Smith  and  Hay  (1961)  laid  the  foundation  for  considering 


cloud  spread  as  a  function  of  turbulent  energy  in  a  "sliding” 
spectral  window.  They  found  both  experimentally  and 
theoretically  that  cloud  spread  will  approximately  follow  a 
linear  growth  relationship 

V-°.22  ix  ,  (2) 


where  i  is  the  total  turbulence  intensity,  and 

2  2.x  ‘/2 

( <u  +v  +w  >) 

1  U 

where  x  is  downwind  distance, 

u,v,w  are  fluctuating  wind  components,  and 
U  is  mean  windspeed. 


(3) 


Mikkelsen  and  Eckman  (198*0  have  recently  given  support  to  these 
results  from  overland  experiments  for  surface  releases  in  short 
to  medium  ranges. 

At  long  distances  from  the  source,  correlation  between 
hypothetical  particles  within  the  cloud  approaches  zero,  and  the 
spread  of  the  cloud  will  behave  like  the  asymptotic  single¬ 
particle  solution, 


I 


Large  scale  motions,  however,  are  often  organized  into  coherent 
vortex  structures  (i.e.  cyclones).  As  a  result,  this  limit  is 
rarely  reached,  and  growth  continues  along  a  more  or  less  linear 
asymptote . 

All  of  these  above  arguments  and  the  rather  linear  shape  of 
the  data  scatter  lead  us  to  a  linear  paramterization .  The 
results  are  shown  in  figures  4-6  and  summarized  in  table  3.  Note 
that  there  is  very  little  unstable  data  at  ranges  greater  than 
3000m.  This  is  primarily  due  to  the  fact  that,  under  these 
conditions,  diffusion  is  enhanced  and  tracer  concentrations  are 
greatly  reduced  far  from  the  source.  This  makes  the  plume 
difficult  to  locate,  and  may  bias  measurements  towards  narrow, 
concentrated  profiles.  The  most  striking  feature  is  that  there 
is  no  clear  distinction  between  P-C  classes  B  and  C  ( GMGO  data), 
or  between  classes  D  and  E  (both  data  sets).  Those  pairs  are  so 
closely  matched  that  we  suggest  only  a  two-class  parameterization 
(also  shown  in  table  3);  one  for  all  unstable  conditions  and  a 
second  for  neutral  to  stable  stratification.  While  turbulence 
intensity  is  not  explicitly  used  in  the  paramterization,  this 
two-class  system  is  supported  by  the  calculated  wind  variances 
from  the  NPS  data  shown  in  table  4.  These  wind  data  include  some 
measurements  from  periods  not  coincident  with  tracer  releases, 
and  do  not  include  some  periods  when  tracers  were  released. 
Nonetheless,  the  largest  change  in  turbulence  (when  moving  from 
one  to  an  adjacent  category)  occurs  between  class  D  and  class  C. 
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DOWNWIND  DISTANCE  (■) 


Figure  4.  Sigma-y  relative  vs.  downwind  distance  with  solid  line 
representing  the  recommended  parameterization  of  table  3.  Symbol 
A  represents  one  data  point,  B  represents  two  data  points,  etc. 

Same  data  as  figure  1,  except  long  range  data  is  included.  Data 
is  a)  Pasquill-Gif ford  stability  class  B,  GMGO  data,  and  b)  class 
C,  GMGO  data. 
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Figure  6.  Same  as  figure  4,  except  a)  Pasquill-Gif f ord  stability 
class  D,  NPS  data,  and  b)  class  E,  NPS  data. 


DATA 

PASQUILL-GIFFORD  CLASS 

SLOPE 

//POINTS 

CMGO 

B 

.0*467 

1  32 

GMGO 

C 

.0357 

1  *40 

CMGO 

D 

.0227 

60’ 

GMGO 

E 

.01  78 

206 

NPS 

D 

.01  78 

215 

NPS 

E 

.0228 

72 

GMGO 

B  +  C  ** 

.0*410 

272 

GMGO  +  NPS 

D+E  ** 

.021  1 

1096 

**Last  2  rows  are  suggested  values. 


Table  3.  Relative  diffusion  parameterization  for  overwater, 
surface,  medium  range  (1-12km)  releases.  Slopes  are  least 
squared  linear  regression  to  data.  S igma-y  relative  = 
(Slope)x(downwind  distance). 


Pasquill-Gifford 
Stab i 1  i tv  Class 


Sigma  Wind  Direction  (deg) 
1  min  ave  10  min  ave 


//  Hours 
Of  Data 


Figure  7.  Wind  direction  sigma  (sigma  theta)  vs.  z/L 
(z  is  the  height,  lDm,  and  L  is  the  Mon in- 
Obukhov  length).  X's  are  1  min  averages, 
and  O's  are  1  hr  averages.  Data  is  from 
NPS.  (see  Schacher,  et  al. ,  1932) 
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4.  MEANDER  (SINGLE  PARTICLE)  PARAMETERIZATION 


For  purposes  of  describing  meander,  information  on  the  motion 
of  the  plume  is  available  only  from  the  NP?  data  set.  Meander  was 
measured  by  calculating  the  distance  between  the  center  of  mass 
for  a  given  profile  and  the  axis  defined  by  the  mean  wind  vector. 
The  method  used  to  calculate  the  mean  wind  vector  was  crucial. 
After  experimenting  with  various  techniques,  we  decided  to  let 
the  average  travel  time  of  the  plume  define  the  proper  averaging 
period  for  the  vector;  this  travel  time  was  roughly  one-half  hour. 
Since  the  wind  was  measured  only  at  the  release  point,  tracer 
profiles  were  correlated  with  average  wind  vectors  in  half-hour 
bins  offset  by  15  minutes  to  account  for  plume  transport. 

The  NPS  experiment  was  designed  so  that  the  plume  was  sampled 
at  nearly  discrete  ranges  f rom  the  release  point.  This  sampling 
procedure  allowed  us  to  segregate  data  into  range  bins.  The 
collection  of  off  axis  center  of  mass  positions  for  a  given  range 
bin  represents  a  probability  distribution  for  the  meander  of  a 
plume  or  puff.  Each  experimental  day  produced  several  distribu¬ 
tions  at  various  downwind  distances  from  the  source.  The  distri¬ 
butions  were  then  analyzed  for  normality,  and  the  standard  devia¬ 
tions  (which  represent  sigma-y  due  to  meander)  were  determined. 

The  data  were  examined  for  dependence  on  windspeed, 
stability  class,  inversion  height,  experimental  day,  and  time  of 
day.  The  day  to  day  variability  was  so  large  that  it  tended  to 
swamp  other  dependencies.  Figures  8-21  show  day  by  day  plots  of 
the  center  of  mass  distributions  interlaced  with  the  respective 
wind  direction  time  histories  and  spectral  plots. 
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TIME  OF  OAY 

Figure  9.  a)  Wind  direction  tine  series,  and  b)  pcwer  spectra 
for  tracer  release  depicted  in  figure  8.  Time  series  data  are 
2  minute  averages  of  wind  direction  sanpled  at  1  Hz.  Spectral 
points  are  smoothed  with  a  5%  sliding  a\erage.  S(n)  is  wind  dir 
ection  pcMer  spectral  density  at  frequency  n.  D  is  the  total 
sanpling  time.  <DIR  >  is  the  mean  wind  direction. 
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Figure  11.  a)  Wind  direction  tine  series,  and  b)  power  spectra 
for  tracer  release  depicted  in  figure  10. 
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Figure  13.  a)  Wind  direction  time  series,  and  b)  power  spectra 
for  tracer  release  depicted  in  figure  12. 
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Individual  data  points  are  marked  with  integer  values  of  the 
measurement  hour  for  a  given  experimental  day.  Note  that  the 
mean  position  of  the  cluster  often  moves  off  the  centerline,  but 
data  are  distributed  evenly  over  the  course  of  a  day.  Large 
variations  in  the  meander  "envelope”  can  be  seen  from  day  to  day, 
and  the  wind  direction  time  series  differ  radically.  The  power 
spectra  also  vary  greatly  from  day  to  day.  Spectral  parameters 
are  summarized  in  table  5.  Spectral  gaps  (defined  in  the  table  5 
caption)  exist  in  most  situations,  with  slopes  of  approximately 
-2  in  the  low  frequency  part  of  the  spectra. 

There  appears  to  be  no  clear  relationship  between  the  total 
energy  in  the  turbulent  and  low  frequency  parts  of  the  spectra. 
Since  meander  is  driven  by  the  low  frequency  energy,  this 
discourages  attempts  at  parameter izat ions  of  meander  based  on 
surface  layer  scaling. 

As  stated  above,  center  of  mass  distributions  were  grouped 
by  experimental  day  and  range  bin.  Variations  of  the  position  of 
the  measurement  aircraft  dictated  range  bins  to  be  2  km  long. 
Relevant  statistics  are  given  in  table  6.  Examination  of  the 
Shapiro-Wilk  statistics  (Shapiro  and  Wilk,  1965)  indicates  that 
the  crosswind  distributions  are  represented  fairly  well  by  normal 
distributions.  This  suggests  that  Gaussian  models  will  be  valid 
even  for  dispersion  caused  by  meander,  given  the  proper  sigma-y. 
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1  Classical  peak  is  the  inverse  of  the  frequency  associated  with 
the  "turbulent  energy  producing"  subrange  maximum  at  the  low 
frequency  end  of  the  -5/3  slope  (inertial  subrange)  region. 

2  Gap  is  the  inverse  of  the  frequency  where  a  minimum  occurs  on 
the  low  frequency  side  of  the  "turbulent  energy  producing" 
subrange . 

3  Slope  is  the  log  (power)  vs.  log  (frequency)  relationship  in 
only  the  "low  frequency"  part  of  the  spectrum.  Gap 
frequencies  are  not  included. 

Table  5.  Summary  of  wind  direction  spectra  for  NPS  data  depicted 

in  figures  9b,  11b,  13b,  15b,  17b,  19b,  and  21b. 
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DATE 

RANGE(km) 

N 

meant  ( m ) 

STD  DEV p ( 

ml  SKEW* 

KURTOSIS  a  We; 

21 Jun82 

2 

25 

-293.0 

149.7 

-0.2861 

-0.3701 

0.9705 

21 Jun82 

4 

26 

-356.5 

255.6 

+0.6759 

+0.0095 

0 . 9554 

22 Jun82 

4 

43 

-1733.7 

392.1 

-0.0101 

-0.6845 

0.9721 

24 Jun82 

2 

1  3 

-80.5 

110.3 

+0.8840 

+0.6898 

0.9077 

24 Jun82 

4 

34 

-678.0 

391  .5 

-0.4567 

+1 . 0885 

0.9590 

25 Jun82 

2 

1  4 

-201  .  3 

50.4 

+0.7334 

-0.8242 

0.8912 

25 Jun82 

4 

38 

-740.9 

166.3 

+0.7936 

+0. 8523 

0.9474 

27 Jun82 

2 

9 

-185.9 

103.3 

-0.0621 

-0.0441 

0.9567 

27  Jun82 

4 

1  7 

-895. 4 

198.8 

-0.1541 

-0  .  1  224 

0 . 967- 

27  Jun82 

6 

4 

-1516.6 

415.4 

+1 .2242 

-0.7663 

'  **  ^  '-4  • 

28 Jun82 

2 

6 

-417.3 

83.4 

+  0.8019 

-0 .1835 

0  .  i  1  ,1 '  j 

28 Jun82 

4 

1  2 

-795.0 

260.9 

-0. 1 653 

+1 .4006 

*•.  n  ■ 

■  -  .  *  M  /  ■  J 

28 Jun82 

6 

6 

-970.6 

160.5 

-1 . 3268 

+1 .7325 

0  .  ••  -  0  - 

28 Jun82 

8 

5 

-1 1 65.7 

467.0 

+2. 1 248 

+4.6043 

■* 

29 Jun62 

2 

3 

-613.9 

477 . 9 

+0.9695 

C  .  ■*  t  1  * 

29 Jun82 

4 

1  2 

-913.0 

1008.9 

+0.6763 

+0 . 0379 

29 Jun82 

6 

9 

-789.9 

966.7 

+0.8435 

-1 .0370 

,  ■■  * 

29 Jun82 

8 

1  1 

-1177.6 

1380.2 

+  1  .2122 

+0 . 7668 

^  i  4 

1  MEAN:  Mean  position  of  plume  center  (as  measured  from  the 

mean  wind  centerline)  of  mass  distribution  for  a 
given  range  bin. 

2  STD  DEV:  Standard  deviation  of  plume  center  of  mass  for  a 

given  range  bin. 

3  SKEWNESS:  Measures  the  tendency  of  deviations  to  be  larger  in 

one  direction  than  the  other.  Values  are  unbounded. 
Positive  values  indicate  a  "tail"  of  values  exists 
in  the  values  larger  than  the  mean,  negative  values 
indicate  a  "tail"  in  smaller  values. 

it  KURTOSIS:  Measures  the  heaviness  of  the  distribution  "tails" 
(-2.0  <  Kurtosis  <  00 )  .  Relatively  large  values 
indicate  normality  assumption  should  be  questioned, 
but  this  statistic  is  unstable  for  a  small  number 
of  samples  which  is  characteristic  of  this 
exper iment . 

5  W :  Shapiro-Wilk  statistic  for  the  center  of  mass 

distribution.  Measures  the  normality  of  the 
distribution  (0.0  <  W  <  1.0).  Values  close  to  1.0 
indicate  more  normal  distributions. 

Table  o.  Meander  statistics  for  NPS  meander  data  .  N  is  the 
number  of  profiles  for  a  given  crosswind  center  of  mass 
distribution. 
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While  we  were  not  able  to  parameterize  meander  with  mean 
surface  layer  meteorological  variables,  direct  measurements  of 


lateral  turbulence  intensity  did  correlate  well  with  the  observed 
meander  standard  deviations.  Taylor  (1921)  predicts  that  when 
the  length  scale  of  turbulence  is  much  larger  than  the  scale  of 
the  plume  (sigma-y  meander  in  this  case)  the  following  simple 
formula  will  apply: 


where  0ym  is  the  single-particle  lateral  standard  deviation,  or 
sigma-y  meander, 

iv  is  the  lateral  turbulence  intensity, 
x  is  downwind  distance. 

This  is  often  referred  to  as  the  "near-field"  approximation  to 
single-particle  diffusion.  For  application  to  the  present 
analysis,  we  define  the  lateral  turbulence  intensity  as 


2.  l/2 


<v  > 

”"u 


(6) 


where  <v^>  is  the  variance  of  the  crosswind  component  measured 

over  the  travel  time  of  the  tracer  (approximately  1/2 
hour)  and  averaged  over  the  entire  day. 

U  is  the  mean  wind  speed 


Figures  22-28  show  equation  (5)  plotted  for  each 
experimental  day  with  the  standard  deviations  of  table  5  (sigma-y 
meander)  depicted  as  data  points.  Also  shown  is  a  least-squared 
linear  fit  to  the  data. 

Note  the  good  agreement  between  this  simple  model  and  the 
observed  data  in  almost  every  case.  The  key  to  success  is 
matching  the  travel  time  of  the  plume  to  the  RMS  window  in  the 
turbulence  intensity  measurement.  Also,  since  the  measured 
variance  is  highly  non-stationary ,  it  is  necessary  to  average 
this  quantity  over  a  full  day. 

The  agreement  between  the  data  and  the  near-field 
approximation  is  also  supported  by  our  spectral  results. 

The  spectra  continue  to  increase  in  power  through  the  entire 
spectral  window  (4-8  hr  )  which  suggests  that  turbulent  motions 
operating  at  these  very  low  frequencies  have  very  long  time 
scales  as  equation  (5)  assumes. 
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Figure  22.  Sigma-y  meander  vs.  downwind  distance.  Error  bars 
are  the  standard  deviations  of  downwind  distance  positions  for 
a  given  range  bin.  Upper  number  is  Shapiro-Wilk  statistic  and 
lower  is  number  of  data  points.  Solid  line  is  best  fit  to  data 
Dashed  line  is  equation  5.  Data  is  NPS,  21  June  1982. 
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Figure  23.  Same  as  figure  22,  except  NPS  data,  22  June  1982. 
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Figure  24.  Same  as  figure  22,  except  NPS  data,  24  June  1982. 
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Figure  25.  Same  as  figure  22,  except  NPS  data,  25  June  1982. 
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Figure  26.  Same  as  figure  22,  except  NPS  data,  27  June  1982. 


5.  CHEMICAL  WEAPONS  HAZARD  FORECAST  PROGRAM  MODIFICATIONS 


Part  of  the  current  research  work  is  to  modify  the  Chemical 
Weapons  Hazard  Forecast  Program  (CWHFP).  This  section  assumes 
the  reader  is  familiar  with  that  program  and  understands  the 
program's  design,  inputs,  and  operation.  EPG  modified  the  24 
May  1984  "test  and  evaluation"  version  of  this  program  supplied 
to  EPG  by  J.  Branum  (1984).  Later  versions  of  this  model  may 
exist,  and  the  following  modifications  can  be  easily  transplanted 
into  such  hybrids. 

The  proposed  modifications  specifically  make  use  of  the 
parameterizations  described  above  and  are  designed  so  as  to 
supply  the  user  with  a  better  understanding  of  the  character  of 
atmospheric  diffusion  processes,  and  hopefully,  more  information 
on  how  to  avoid  hazards.  The  previous  model  used  meteorological 
inputs  to  calculate  a  stability  category,  and  then  calculate  one 
hour  averaged  plume  dimensions  based  on  stability  dependent 
sigma-y  and  sigma-z  algorithms.  This  version  uses  the  same 
stability  scheme,  but  replaces  the  one-hour  sigma-y  functions 
with  the  relative  diffusion  parameterization  presented  above. 
Sigma-z  parameterizations  are  unchanged.  The  resultant 
isopleths  represent  either  "dosage"  values  from  an  instantaneous 
"burst"  release  or  a  concentration  "snapshot"  from  a  continuous 
plume . 
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In  order  to  present  the  meander  of  a  plume  (offaxis 
excursions  of  a  puff)  the  single-particle  parameterization 
developed  above  is  used  to  obtain  a  "meander  envelope."  This 
envelope  is  simply  the  2  sigma-y  position  of  the  plume  or  puff 
center  of  mass  distribution  representing  a  95$  probability  of 
impact.  This  envelope  is  superimposed  on  the  instantaneous 
isopleths  in  order  to  show  the  combined  impact  of  the  two 
diffusion  processes.  An  artificial  "ripple"  is  convolved  into 
the  instantaneous  solution  to  offer  the  user  a  visualization  of 
meander.  The  frequency  and  amplitude  of  this  ripple  has  no 
physical  basis,  and  is  only  installed  to  help  the  user  understand 
the  nature  of  the  meander  process. 

Figures  29-36  show  several  examples  of  model  output.  Labels 
display  lethal  dose  levels  on  the  instantaneous  isopleths  and  a 
footnote  explains  the  meander  envelope.  Instruction  manuals 
should  explain  the  arbitrary  nature  of  the  ripple  in  order  to 
avoid  false  conclusions  based  on  these  plots.  Note  the  change  in 
the  frequency  of  the  ripple  as  the  meander  envelope  shrinks  (see 
figures  29,  30,  31  or  32,  33,  3^).  This  change  again  was 
arbitrary  ,  and  was  only  intended  to  suggest  the  possibility  that 
the  important  eddy  sizes  decrease  as  the  meander  decreases.  The 
amplitude  of  the  ripple,  again,  was  arbitrary  (selected  to  be  1 
sigma-y  meander)  and  suggests  that  the  instantaneous  position  of 
the  cloud  is  somewhat  random  and  rarely  reaches  the  meander 
envelope  bounds.  Also  note  the  change  in  the  instantaneous  cloud 
dimensions  with  a  change  in  stability  class  (figures  29,  30,  31, 
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Figure  29.  Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  unstable  atmosphere,  large  meander. 
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Figure  30.  Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  unstable  atmosphere,  moderate  meander. 
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Figure  31.  Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  unstable  atmosphere,  small  meander. 
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Figure  32.  Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  stable  atmosphere,  large  meander. 
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Figure  33.  Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  stable  atmosphere,  moderate  meander. 


CvC*»  V*r*r 


CHEMICAL  UEPPON  HRZfiPD  f  0RECA9T  PROGRAM  -  SNAP  XX.  >' 


RINGS  -  4  00  YARDS  BERPENGS  IN  DEG  TRUE  PLOT  FORMAT  ♦! 


TERRAIN  TYPE  OPEN-SEA 

MEAN  MIND  3  KTS  FROM  273  DEC  TRUE 

SIGMA  WIND  DIRECTION  12  DEC 

MEAN  RIP  TEMP  10  DEC  C 

MEAN  SEA-SURFACE  TEMP  17  DEG  C 

RELATIVE  HUMIDITY  70  H 

SURFACE  MIXING  LRYEP  HT  0  METERS  BASED  ON  #•  DEFRULT  •* 

STABILITY  CRTRGORY  B  MODIFIED  PRSQUILL 

MUNITION  TYPE  MK  1  16-SIZE  BOMB'M  I SS  UE 

source  type  point-burst 

SOURCE  SIZE  i«ff«ctts/*>  88  KG 

SOURCE  RRTE  INSTANTANEOUS 


contour  LABEL  POTENTIAL  CASUALTY  effects  approx  MAX 

1  DOSE  -AGENT  y  _ •  Wl  ThQUT  PROTECT  ION) _  RANGE 

-  LD30-GD  30*.  DEATHS  -  MOST  INCAPACITATED  1301  YARDS 

-  L  D 1 -CD  r.  DEATHS  -  MANY  INCAPACITATED  2872  YARDS 

-  LD0-GD  NO  DEATHS  -  SOME  INCAPACITATED  3846  YARDS 

OUTSIDE  -THICK'  CONTOUR  IS  S3*.  PROBABILITY  THAT  CENTER  OF  MASS  HILL  CROSS 


Figure  35.  Chemical  Weapons  Hazard  Forecast  Program  output 
sample  —  light  wind  conditions. 


CHEMICAL  WEAPON  HAZPiPD  roRECPST  PROGRAM  -  SNAP  XX.  X 


PINGS  -  1000  yhRDS 


TERRAIN  T  t P£ 

HERN  MIND 

SIGMA  WIND  DIRECTION 
MENU  AIR  TEMP 
MERN  SEA-iURFACE  TEMP 
RELATIVE  HUMIDITY 
SURFACE  MIXING  LAYEP  HT 
STABILITY  catagory 
MUNITION  TYPE 
SOURCE  T.PE 

SOURCE  SIZE  *vt ) 

SOURCE  RATE 


BERRENGS  IN  DEG  TRUE 


OPEN-SER 
20  KTS 


PLOT  FORMAT  #1 


73  DEG  TRUE 


0  METERS  BASED  ON  DEFAULT  •• 

D  MODIFIED  PASAUILL 

MM16-SIZE  BuMB'MISSILE 

POINT-BURST 

86  KC 

INSTANTANEOUS 


vs.  32,  33,  34)  or  wind  speed  (figures  35  vs.  36)  as  would  be 
qualitatively  predicted  by  a  Gaussian  model.  This  model  requires 
one  additional  user  input  over  those  needed  in  the  previous 
version;  the  standard  deviation  of  the  wind  direction  used  in  the 
meander  calculations.  As  stated  in  the  previous  sections,  the 
averaging  time  must  be  the  cloud  travel  time  to  the  distance  of 
interest.  Also,  the  value  used  should  be  the  average  of  several 
such  standard  deviations  obtained  over  a  significant  fraction  of 
the  day.  The  sampling  time  is  not  critical  ,  1.0  -  0.1  Hz  would 
be  sufficient. 

Code  modifications  are  listed  in  figures  37-40.  Figure  37 
shows  both  the  turbulence  intensity  input  and  the  revised  sigma 
parameter  table  (relative  diffusion).  Figure  38  shows  where  the 
meander  envelope  is  calculated  while  figures  39-40  list  the 
plotting  routine. 


Mean_w i nd_sp_k  t *5 
Mean_w  1  nd_d i r  =  275 
! 

Mix  1  aver  ht  *0 


!  KTS 

!  DEG  TRUE 
!  METERS 


!  *****  NEW  1985  ********* 

Turb_i  ntensM2*PIM80  1  turbulence  intensity  (raflian:  )  input 

i  ************************ 

i 

DIM  Ml_ht_est_meth*C30I 

Ml  ht  est  met h** " **  DEFAULT  **" 


3790  !  LOAD  ARRAY  CONTAINING  VALUES  OF  Ap->Dp 

3800  ' 

3810  FOR  POM  TO  6  !  WHERE  P0  »  NUMERICAL  REFERENCE  TO  STABILITY  CLASSES 

3820  i  <  1 «  "  A  "  ->  6*  " F "  ) 

3830  READ  Mat r i x_ap< P0 ) , Mat r i x  bp ( P0 ) , Mat r i x  c p ( P0 > , Mat r i x  dp(P0> 

3840  NEXT  P0 
3850  ' 

3860  '  *******  NEW  1985  ********* 

3870  i 

3880  DATA  0,0, 0,0  !  --  CAT  “A"  (NOT  USED  OVER  THE  OCEAN) 

3890  DATA  0.0410,1.00,0.32,0.75  !  --  CAT  “B" 

3900  DATA  0.0410,1.00,0.32,0.70  •  --  CAT  "C"  sigma*  are  for  relative  diffusion 

3910  DATA  0.0211,1.00,0,16,0.65  1  --  CAT  "D"  in  lateral  direction  only 

3920  DATA  0.0211,1.00,0.10,0.62  !  --  CAT  "E" 

3930  DATA  0,0, 0,0  !  --  CAT  "F-V'CNOT  DEFINED  OR  USED  BY  THIS  MODEL) 

3940  ' 

3950  !  ************************** 

3960  ' 


5710  !  ****************************#****************************«************* 
5720  i 

5730  Pt_«*por_nomx:  '  CALC  MODEL  FOR  POINT  SOURCES  OF  GASEOUS  AGENTS, 

5740  !  AND  NO  REFLECTIVE  MIXING  LAYER  CAP 

5750  ! 

5760  ! 

5770  ! 

5780  <  *****  NEW  1985  **** 

5790  DIM  Meand_x< 100) ,Meand_y< 180) 

5880  DIM  Bat a_set_x < 3 , 1 00 ) , Dat a_se t_y < 3 , 1 00 >  1  10O  point*  ms.  20  pts 

5810  Points_per  set=100  !  in  old  version 

5820  !  ******************** 

5830  ! 

5840  M_to_yds=39. 37 '36  ‘  CONVERSION  FACTOR  USED  IN  EQUATIONS  BELOW,  DEFINED 

5850  !  AGAIN  HERE  TO  PERMIT  TESTING  OF  PARTIAL  PROGRAM 

5860  • 

5870  FOR  K=1  TO  Num_dat a_set s 

5880  Part  i  a1_l  < K ) *Sour c e_s i  ze (. K > / < P  I  *60*Mean_w i  nd_sp_m*Dose_val  (K)  > 

5890  Max_r ange_x  <K >  =  < Part i al _1 <K >x < Ap( K) *Cp<K >  >  > ~< 1 x< Bp<K)+Dp<K> ) ) 

5900  Dat a_set_y cK, 0)*0 

5910  X_i ncrement “Max_range_x<K >xPo i nt s  per_set 

5920  IF  X_i nc r ement  =0  THEN-X_i nc rement  =  1  !  ADDED  TO  PREVENT  BLOWUP  ON  BHD  DATA 

5930  ! 

5940  FOR  3=1  TO  Po l nt s_per_*et 
5950  X  = J*X_i ncrtment 
5960  Data  set_x<K, J>=X 

5970  S l gma_yx  =  Ap  <  K)-*X^BpOO  !  ASSUMES  NO  INITIAL  CLOUD  SI2E  FOR  CONSERVATISM 
5980  Si  gma_zx  =  Cp<K)*X'Dp(.K) 

5990  Partial  2<K )=2*LOG<Part i al  lOOxCSigma  yx#Sig#ia  zx>> 

6800  IF  Partial  2<K><0  THEN  Partial  2<K)=0  !  PREVENT- POSS I BLE  ERROR  IN  LATER  SQR 

CALCULATIONS 

6010  Data  set  yCK, 3)=Si gma  yx*SQR<Part i al  2<K>) 

6020  NEXT- J 
6030  NEXT  K 
6040  ! 

6050  !  *****  NEW  1985  ********* 

6068  RAD 

6070  FOR  J=1  TO  To i nt s_per_set 
6080  X=3*X_i ncrement 

6090  IF  X>Max_range_x<Num_data_sets)*C0S<2*Turb_intens)  THEN  GOTO  Close 
6100  Meand_y<  3 ) =X*T AN < 2*Turb_i nt ens )  !  meander  envelope  is  based 

6110  Meand_x<3)=X  !  only  on  turbulence  intensity  and 

6120  GOTO  310O  !  selected  to  be  the  2  sigma  value  (95H) 

6130  ! 

6140  Close:  Meand_y< 3 ) *SQR < Max_r ange_x < Num_dat a_se t s > ~2-X~2 )  !  stop  envelope 
6150  Meand_x< J)*X  !  at  maximum  range  and  close  with  partial  circle 

6160  3100:  NEXT  J 

6170  MAT  Meand_x=Meand_x*<M_to_yds) 

6180  MAT  Meand_y=Meand  y*(M  to  yds) 

6190  ! 

6200  '  *********************** 


Figure  38.  Abbreviated  listing  of  CWHFP  showing  calculation 
of  the  meander  envelope.  Points  per  data  set  may  be  adjusted 
for  speed  or  better  resolution. 


ei9o  1  *****  new  1965  immmt 

8 2 O 0  DIM  Seal  ed_»  (100>  ,  Seal  ed_y(  100)  1  LIMITS  EftSED  ON  "POINTS_PER_SET" 

8218  DIM  Corit  our_l  ibl_x\3),Contour  1abl_y(3)  !  LIMIT  BASED  ON  NUM_DRTA_SETS 
8 2 2  0  1 

8230  FOP  Point  num=0  TO  Points  per  set  1  PO I NTS_PEP_SET  preset  above 
8240  Sc  al  ed_*tPoi  nt  _num)=Meand_xv  Poi  nt_num > *Sc al  e_f  ac  t  or  !  scale  meander 
8250  8c al e d_y> Po i nt _num > =Meand_y (Po i nt_num) *Sc al e_f ac t or  !  envelope 
6260  NEXT  Point _num 
6270'  GOSUB  Plot_meander 
6280  • 

8290  >  ************************** 

8300  1 

8310  FOP  Set_num=l  TO  Num_dat a_set s  1  NUM_DATR_SETS  PRESET  TO  “ 3 “ ,  ABOVE 

8  320  1 

8330  FOP  Point  num  =  0  TO  Po i nt s_per_set  1  PO I NTS_PER_SET  PRESET  ABOVE 
8  340  Scaled  x  (  Po  i  nt  _num  >  *Dat  a_se  t  ._x(  Set  _num ,  Po  i  nt _num  >  *Sc  al  #_f  actor 
8  350  Seal  ed_u  (Point _num>=Dat  a_set_y ( Set  _num,  Po  i  nt_nui»  > *Sc  al e_f  ac  t  or 
8360  NEXT  Point_nuni 
8  370  GOSUB  F'l  ot  _contour 

8  380  C  ont  our _1 ab 1  x  < Se t _num >  =  Sc  al ed_x (Point s_per_set  > 

8  390  1 

8400  '  tittstttt  new  1985 

8410  RhD 

8420  P= D at a_se t _x < Se t  num , Po i nt s  per  se t > xDat a_set_x ( Num_dat a_set s , Po i nt s_per_s 
e  '  ■  1  d i mens i on  1  ess  distance  from  source  [0,11 

8430  Ome ga=48*P I ' 1 80  'Turb_ i nt ens  !  ARBITRARY  frequency  of  ripple  oscillation 
8440  A=S: al ed_x(Po i nt s_per_set >*TAN< Turb_i nt ens)  1  ARBITRARY  ripple  amplitude 
8450  R i pp 1 e=A*S I N( Omega*P I *R >  1  offset  for  y  position 

8460  Con  t  our  _1  ab  1  _y  C  Se  t  _num  >  “Sc  al  ed_y  <  Po  i  nt  s_per_se  t  )  +R  i  pp  1  e  !  add  offset  to  la 
be  1  position 
8470  DEG 

8480  *  ** ** ******************************************* 

8490  ' 

8500  NEXT  Set_num 
8510  GOTO  Exit_data  plot 
8520  ' 

8530  '  *******  NEW  1985  ******************************  (NEW  SUBROUTINE) 

8540  ' 

6550  F 1  o’  _n,e  ander :  ’  SUB  to  plot  meander  envelope 

8560  LINE  TYPE  1 

8570  Max_o  "er st r i k e* 1  !  same  as  "old"  version 

8580  X_of f  set  =0 

8590  1  TOP  HALF . 

8600  FOR  St r i k ec ount =0  TO  Max_over st r i k e 
8610  Y_cf f set*. 33«Strikecount 

8620  X_of t s e t = . 33*St r i k ec ount  !  strikecounts  thicken  line 
8630  MOVE  Sc  a  1 ed_x < 0 ) , Sc  a  1 ed_y < 0 ) 

8640  FOR  t  =0  TO  Po i nt s_per_set 

8650  DRAW  Sc  a  1 ed_x < K > +  X_of f set , Sc al ed_y < K )♦ Y_of f se t 
8660  NEXT  K 

8670  1  SCALED_X,Y  now  contain  meander  envelope 

8680  1  PLOT  •'BOTTOM"  HALF 

8690  MOVE  Sc  a  1 ed_x < 0 ) , Sc al ed_y < 0 ) 

6700  FOP  K=0  TO  Po i nt s_per_set 

8710  DPfiw  Scaled_x(K>fX_off*et,-(Sca)*d_y(K>4Y_offs*t> 

6720  NEXT  k 

8730  NEXT  Strikecount 
8740  RETURN 

8750  '  *************************** 

Figure  19.  Abbreviated  listing  of  CWHFP  showing  meander  envelope 
plotting  routine  and  contour  label  locating  scheme. 


8770  !  tttmi  NEW  1985  fttftttttt  (.REVISED  SUBROUTINE) 

8780  Plot  contour:  1  plot  contours  of  the  scaled  data 
8790  RAB~ 

8800  LINE  TYPE  1 

8810  0m« ga“48*P I •' 180 ' Turb  inters  !  ARBITRARY  frequency  of  ripple  oscillation 

8820  > 

8886  1  PLOT  BEGINNING  WITH  "TOP"  HALF  OF  CONTOUR 
8340  !  <  OVERSTRIKE  STUFF  OMITTED  ) 

8850  MOVE  Sc al ed_x t 0 > , Sc al ed_y <0 > 

88G0  FOR  1=0  TO  Po i nt s_per_set 

8870  R=Bat a_se t_x < Set_num , K > ' Bat a_se t _x< Num_dat a_se t s , Po i nt s_per_set )  !  dimensi 

onless  distance  from  source 

8880  A=Sc al ed_x ( K ) *TAN< Turb_i nt ens >  !  ARBITRARY  amplitude  of  ripple  selected  to 
be  1  sigma 

8890  Ri pp 1 e=A*S IN(0mega*P I *R >  !  offset  in  y  direction 
8900  DRAW  Scaled  x<K), Scaled  y<K>+Ripple  !  add  offset 
8910  NEXT  K 
8928  1 

8930  !  PLOT  "BOTTOM”  HALF  ....  most  same  as  “top"  half 
8940  • 

8958  MOVE  Scaled_>(0),Scaled_y(0) 

8960  FOR  K  =0  TO  Po i nt s_per_set 

8970  R=Bat  a_set  _x < Set_num, K )  .'Dat a_set_x<Num_dat  assets .Point  s_per_set  > 

8980  A=Scaled_x(K >*TAN< Turb_i nt ens > 

8990  Pi ppl e=A*SIN<Omega*PI*R> 

9000  DRAW  Scaled  x ( K ) , - < Sc al ed  y(K)-Ripple)  1  subtract  offset 
9810  NEXT  K 
9020  DEG 

9030  RETURN  !  END  OF  GOSUB  PLOT  CONTOUR 
9040  !  *t*ft**S*t*S$fSS***S*$*Sf*SS**$SS$St 


Figure  40.  Abbreviated  listing  of  CWHFP  showing  plotting  of 
hazard  contours  based  on  relative  diffusion  parameterization. 
"Wiggles"  are  a  result  of  a  SIN  wave  imposed  on  the  Gaussian 
plume  model  solution. 


6.  SUMMARY  AND  FUTURE  WORK 


The  presented  work  successfully  parameterizes  relative 
lateral  diffusion  from  mean  meteorological  quantities,  verified 
via  two  independent  overwater  tracer  data  sets.  Single- particle 
diffusion  is  successfully  parameterized  only  through  direct 
measurements  of  the  lateral  turbulence  intensity.  Once  this 
quantity  is  known,  however,  the  parameterization  is  well  behaved. 

Future  work  should  attempt  to  correlate  this  lateral 
turbulence  with  large  scale  synoptic  features,  mesoscale 
phenomena,  radiosondes  profiles,  or  other  more  easily  measured 
quantities.  The  GMGO  researchers  are  presently  being  petitioned 
for  positioning  information  on  their  experiments  so  that  a  more 
thorough  verification  of  the  single-particle  parameterization  can 
proceed  in  the  future. 

The  model  is  presently  designed  for  only  medium  ranges,  but 
could  tie  extended  to  greater  distances.  This  extension  would 
make  results  highly  dependent  on  inversion  height.  Some 
prognostic  or  diagnostic  estimate  of  this  quantity  could  be 
incorporated  into  the  model.  Longer  range  overwater  experiments 
.'■ho  i  1  !  he  researched  for  plume  pa ramter i za t i ons  at  these 
1  i  stunner. . 
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ATTN:  LIBRARY,  COOE  2620 
WASHINGTON,  DC  20390 


COMMANDING  OFFICER 
NAVAL  OCEAN  RSCH  &  DEV  ACT 
NSTL,  MS  39529-5004 


COMMANDING  OFFICER 
NAVAL  OCEANOGRAPHIC  OFFICE 
BAY  ST.  LOUIS 
NSTL,  MS  39522-5001 


CHAIRMAN 

OCEANOGRAPHY  DEPT. 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 


NAVAL  POSTGRADUATE  SCHOOL 
METEOROLOGY  DEPT. 
MONTEREY,  CA  93943-5000 


LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 
MONTEREY,  CA  93943-5002 


COMMANDER 

NAVAIRSYSCOM  (AIR-330) 
WASHINGTON,  DC  20361-0001 


CHIEF,  ENV.  SVCS.  DIV. 
OJCS  ( J-33) 

RM.  2877K,  THE  PENTAGON 
WASHINGTON,  DC  20301 


COMMANDING  OFFICER 
OFFICE  OF  NAVAL  RESEARCH 
1030  E.  GREEN  ST. 
PASADENA,  CA  91101 


COMMANDER 

OCEANOGRAPHIC  SYSTEMS  PACIFIC 
BOX  1390 

PEARL  HARBOR,  HI  96860 


COMMANDING  OFFICER 
FLENUMOCEANCEN 
MONTEREY,  CA  93943-5005 


DIRECTOR  OF  RESEARCH 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 


NAVAL  POSTGRADUATE  SCHOOL 
OCEANOGRAPHY  DEPT. 
MONTEREY,  CA  93943-5000 


PRESIDENT 

NAVAL  WAR  COLLEGE 

GEOPHYS.  OFFICER,  NAVOPS  DEPT. 

NEWPORT,  RI  02841 


COMMANDER 

NAVAIRSYSCOM 

MET.  SYS.  DIV.  (AIR-553) 

WASHINGTON,  DC  20360 


NAVAL  DEPUTY  TO  THE 
ADMINISTRATOR,  NOAA 
ROOM  200,  PAGE  BLDG.  #1 
3300  WHITEHAVEN  ST.  NW 
WASHINGTON,  DC  20235 


OFFICE  OF  NAVAL  RESEARCH 
SCRIPPS  INSTITUTION  OF 
OCEANOGRAPHY 
LA  JOLLA,  CA  92037 


COMMANDER 

NAVAL  OCEANOGRAPHY  COMMAND 
NSTL,  MS  39529-5000 


SUPERINTENDENT 
LIBRARY  REPORTS 
U.S.  NAVAL  ACADEMY 
ANNAPOLIS,  MD  21402 


SUPERINTENDENT 

NAVPGSCOL 

MONTEREY,  CA  93943-5000 


NAVAL  POSTGRADUATE  SCHOOL  (75) 
PHYSICS  &  CHEMISTRY  DEPT. 
MONTEREY,  CA  93943-5000 


COMMANDER  (2) 

NAVAIRSYSCOM 
ATTN:  LIBRARY  (AIR-723D) 
WASHINGTON,  DC  20361-0001 


COMMANDER,  SPACE  &  NAVAL 
WARFARE  SYSTEMS  COMMAND 
ATTN:  CAPT.  K.  VAN  SICKLE 
CODE  06G,  NAVY  DEPT. 
WASHINGTON,  DC  20363-5100 


COMMANDER 

NAVAL  SEA  SYSTEMS  COMMAND 
ATTN:  LT  S.  PAINTER 
PMS-405-24 

WASHINGTON,  DC  20362-5101 


COMMANDER 

NAVAL  SHIP  RSCH  &  DEV.  CENTER 
CODE  5220 

BETHESDA,  MO  20084 


COMMANDER 

NAVAL  SURFACE  WEAPONS  CENTER 
ATTN:  CODE  44 
DAHLGREN,  VA  22448-5000 


USAFETAC/TS 
SCOn  AFB,  IL  62225 


COMMANDER  &  DIRECTOR 
AHN:  DELAS-D 
U.S.  ARMY  ATMOS.  SCI.  LAB 
WHITE  SAND  MISSILE  RANGE 
WHITE  SANDS,  NM  88002 


COMMANDER  &  DIRECTOR 
ATTN:  DELAS-AS 
U.S.  ARMY  ATMOS.  SCI.  LAB 
WHITE  SANDS  MISSILE  RANGE, 
NEW  MEXICO  88002 


DIRECTOR  (12) 

DEFENSE  TECH.  INFORMATION 
CENTER,  CAMERON  STATION 
ALEXANDRIA,  VA  22314 


HEAD,  ATMOS.  SCIENCES  DIV. 
NATIONAL  SCIENCE  FOUNDATION 
1800  G  STREET,  NW 
WASHINGTON,  DC  20550 


ROGER  HELVEY,  CODE  3253  (4) 

PACMISTESTCEN 

PT.  MUGU,  CA  93042 


COMMANDER 

NAVOCEANSYSCEN 

DR.  J.  RICHTER,  CODE  54 

SAN  DIEGO,  CA  92152-5000 


COMMANDER 

NAVSURFWEACEN,  CODE  R42 
DR.  B.  KATZ,  WHITE  OAKS  LAB 
SILVER  SPRING,  MD  20910 


COMMANDING  OFFICER 
AMPHIBIOUS  WARFARE  LIBRARY 
NAVAL  AMPHIBIOUS  SCHOOL 
NORFOLK,  VA  23521 


COMMANDER 

ARMY  DEFENSE  &  READINESS  COM 
AHN:  DRCLDC 
5001  EISENHOWER  AVE. 
ALEXANDRIA,  VA  22304 


COMMANDING  OFFICER 
U.S.  ARMY  RESEARCH  OFFICE 
AHN:  GEOPHYSICS  DIV. 

P.O.  BOX  12211 
RESEARCH  TRIANGLE  PARK,  NC 
27709 

COMMANDER  &  DIRECTOR 
U.S.  ARMY  ATMOS.  SCI.  LAB. 
AHN:  DELAS-AF 
WSMR,  NEW  MEXICO  88002 


DIRECTOR,  TECH.  INFORMATION 
DEFENSE  ADV.  RSCH  PROJECTS 
1400  WILSON  BLVD. 

ARLINGTON,  VA  22209 


DIRECTOR 

FEDERAL  EMERGENCY  MANAGEMENT 
AGENCY  (FEMA) 

WASHINGTON,  DC  20472 


STUART  GATHMAN,  CODE  4117 
NAVAL  RESEARCH  LAB 
WASHINGTON,  DC  20375 


COMMANDER 

NAVAL  WEAPONS  CENTER 
DR.  A.  SHLANTA,  CODE  3331 
CHINA  LAKE,  CA  93555-6001 


DIRECTOR 

NAVSURFWEACEN,  WHITE  OAKS 
NAVY  SCIENCE  ASSIST.  PROGRAM 
SILVER  SPRING,  MD  20910 


COMMANDER 

AWS/DN 

SCOn  AFB,  IL  62225 


DEPT.  OF  THE  ARMY 
CHIEF  OF  ENGINEERS  OFFICE 
ATTN:  DAEN-RDM 
WASHINGTON,  DC  20314 


DIRECTOR 

LIBRARY,  TECH.  INFO.  CEN. 
ARMY  ENG.  WATERWAYS  STN. 
VICKSBURG,  MS  39180 


COMMANDER/DIRECTOR 

US  ARMY  ATMOS.  SCIENCE  LAB. 

AHN:  DELAS-AT-0 

WHITE  SANDS  MISSILE  RANGE,  NM 

88002 


COMMANDANT 

U.S.  COAST  GUARD 

WASHINGTON,  DC  20226 


DR.  MARVIN  DICKERSON 
L-262,  LLNL 
P.O  BOX  808 
LIVERMORE,  CA  94550 


JIM  BRANUM 

2682  OLIVESTONE  WAY 

SAN  JOSE,  CA  95132 
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